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Abstract: Cognitive maps are a vital tool that can be used for knowledge representation and 
reasoning. Fuzzy Cognitive Maps (FCMs) are popular soft computing techniques used to model large 
and complex systems, and they can aid in explainable artificial intelligence (AI). FCMs, however, 
cannot model the indeterminacy that arises in a system due to various uncertainties. Neutrosophic 
Cognitive Maps (NCMs), upgraded FCMs that could model indeterminacy, were introduced to 
address this issue. NCMs are a generalization of FCMs, a field of cognitive science firmly based on 
neural networks. NCMs have been used to solve a wide range of problems. NCMs were introduced 
in 2002, and even after 20 years, NCMs do not have any supportive software, package, toolbox, or 
visualization software like FCMs. The main reason for the absence of dedicated software is due to the 
indeterminacy concept 'I' and how it has to be handled. This paper presents the dedicated Python 
package created for handling the functioning of NCMs. The modelling software presented in this 
paper aids in visualizing the NCMs as a signed digraph with indeterminacy that is a directed signed 
neutrosophic graph. This package implements a sample case study using NCMs. 


Keywords: Neutrosophy; Neutrosophic Cognitive Maps; Python Package; Visualization of NCMs. 


1. Introduction 


Fuzzy theory is a branch of mathematics that deals with vagueness and uncertainty in decision- 
making [1]. Fuzzy sets and logic model complex problems involving imprecise terms or partial truths. 
It has many applications in engineering fields, the healthcare sector, economics, and social science, 
which pertain to real-world problems. Fuzzy logic and its models have many applications in various 
fields, such as engineering, artificial intelligence, medicine, economics, and social problems. It can 
help model multifaceted problems that involve human knowledge, preferences, or emotions [2]. 

A fuzzy cognitive map (FCM) [3] represents a mental landscape within which the connections 
between the nodes (e.g., events, concepts, resources, or attributes) are used to compute the “strength 
of impact" of these elements. FCMs are signed fuzzy digraphs introduced by Bart Kosko [3]. 

FCMs have been used to analyze several socio-economic, healthcare, and decision-making 
problems. The applications and extensions of FCMs are vast and widely researched; a few are 
presented here. 

In [4], the authors explore using FCMs as an agency for collective decision-making and how 
FCMs can capture the cognitive models and group beliefs of different stakeholders. FCMs were used 
as a learning assessment tool in [5] to understand the planning of children by stimulating cognitive 
function. 
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In [6], the Multi-Agent Genetic Algorithm (MAGA) is proposed to optimize convergence error 
for learning FCMs. A GIS-dependent crisis management tool to predict earthquakes in Tehran using 
FCMs was proposed in [7]. 

An evolutionary algorithm, known as IBMTEA-FCM, was proposed in [8] for learning large- 
scale FCMs. A qualitative analytical method using FCMs to specify causal-effect links between the 
interdependent SDGs by considering the long-term effect of COVID-19 was presented in [9]. 

In [10], the PRescriptiVe FCM (PRV-FCM) was introduced, based on FCMs and metaheuristic 
algorithms, to develop prescriptive models. In [11], three federated learning approaches were 
combined with FCMs for mortality prediction and treatment prescription in severe dengue cases. 

Iran’s population’s health was analyzed [12] using FCMs. Since the concept of health is a 
complex and comprehensive system, other sector policies profoundly affect health. Borrero- 
Dominguez and Escobar-Rodriguez [13] proposed a decision support system for crowdfunding using 
FCMs. The various extensions of FCMs have been systemically reviewed by [14]. 

Several software packages are available for modelling FCMs. One such software, Mental 
Modeler, helps build FCMs intuitively and easily. After creating the models, decreasing or increasing 
the model’s elements allows us to examine various change tactics. FCMexpert is a software tool for 
FCM-based scenario analysis and pattern classification presented in [15]. Over 10 FCM extensions 
were handled by supporting interoperability in the FCM extensions in [16]. 

Python packages are also available for modelling FCMs. FCMpy [17] is a recently introduced 
open-source package for building and analyzing FCMs. Notably, FCMPy allows simulating system 
behaviour using qualitative data to create fuzzy causal weights, applying ML algorithms to modify 
the FCMs matrix to aid in classification, and executing scenario examination by simulating theoretical 
interventions. 

The package also helps apply ML algorithms (e.g., nonlinear and active Hebbian learning, 
deterministic learning, and genetic algorithms) to adjust the FCM weight matrix. 

Neutrosophy is a branch of philosophy investigating neutralities’ origin, nature, and scope and 
their interactions. Florentin Smarandache introduced neutrosophy in the 1990s [18]. Neutrosophy 
regards a proposition, hypothesis, concept, event, or entity depending on the modelling. 
Neutrosophy is the basis of the neutrosophic set, logic, probability, and statistics. Indeterminacy ("I") 
is a concept in neutrosophy that measures the degree of neutrality or uncertainty of a proposition, 
event, theory, entity, or concept. 

Neutrosophic Cognitive Maps (NCMs) are an extension of FCMs that can handle indeterminate 
relationships between two concepts, obtaining more significant and sensitive results. It was 
introduced in [19] to analyze diverse social issues. NCMs have been modelled considerably on the AI 
focus to mimic the thinking-humanly approach. Here, it is unsupervised data and has a limited set of 
features. 

Over the past two decades, many investigators have utilized NCMs to analyze diverse problems 
like situation analysis [20], pest analysis [21], transgressions against people experiencing 
homelessness [22], and imaginative play in children [23]. FCMs and NCMs on COVID variants were 
compared in [24]. SWOT analysis and NCMs were combined to analyze organic farming in India [25]. 
Al-Subhi et al. [26] proposed triangular NCMs and used them in multistage decision-making with a 
use case of evaluation. NCMs and cloud data were used in [27] to detect violence, and several datasets 
were used. NCMs and FCMs were compared in this analysis, and it clearly states that NCMs are 
better at handling indeterminacy than FCMs. 

Dynamic NCMs [28], enhanced cuckoo search, and ensemble classifiers were presented for 
acquiring the profile of gene expression and differentiating between the individuals affected by 
rheumatoid arthritis and possible control subjects. Bhutani et al. [29] proposed a technique combining 
pest analysis based on fuzzy and neutrosophic logic to analyze the food industry. 
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In [30], the various factors impacting the paper-packaging industry are analyzed to provide a 
notional representation using NCMs since sustainable supply chains can be attained with repeated 
product use and recycling. 

In [31], NCMs were used to analyze the various causes and effects that lead to violent behaviour. 
NCMs were used to determine the elements that enable proper decision-making to confirm a precise 
diagnosis of conversion disorder [32]. The various factors that affect homeless people were analyzed 
using NCMs in [33]. A neutrosophic sociogram-based NCM approach was introduced in [34]. FCMs 
and NCMs have been used in health care to analyze dengue fever [35]. 

NCMs have been applied in various domains such as health, agriculture, engineering, social 
problems, business, law, environment, and medicine. The substantial advantage of NCMs over other 
cognitive maps is their capacity to capture data realistically and consistently represent expert 
opinions, making them a valuable tool for decision-making strategies where there is an advanced 
degree of indeterminacy or uncertainty. NCMs can be constructed either by a data-driven approach 
or experts’ opinions. 

NCMs have not been integrated with machine-learning algorithms like FCMs. So, little research 
combines various machine learning algorithms that have been utilized in adjusting weights in an 
NCM, like in FCM. Despite the various applications of NCMs, there is no dedicated software or 
Python package for NCMs. This paper presents a dedicated software and Python package for 
constructing, analyzing, and visualizing NCMs. 

The dedicated modules of our modelling software function in the following way: 

i. | Generating the neutrosophic graph and the related connection matrix using expert opinion 
e From linguistic terms. 
e From edge weights. 
e From a file as input from the user. 
ii. Visualizing the NCMs as a neutrosophic digraph. 
iii. | Simulating the NCMs using various state vectors. 
iv. Analysis of various case scenarios for given NCMs. 


The paper is organized as follows: Section 2 recalls the workings of NCMs and their construction. 
Section 3 provides the software-based visualization and simulation of NCM for a case study. The 
conclusions and results are presented in the last section, together with suggestions for future research. 


2. Working of Neutrosophic Cognitive Maps (NCMs) 


When data is unsupervised and the association between two concepts is indeterminate, the 
indeterminacy can be captured by using neutrosophy. [19] introduced the concept of indeterminacy 
in FCMs, called NCMs. The basic properties of NCM are recalled to make this section self-contained. 

NCM is a digraph with concepts as nodes and their causal relationships as edges. These concepts 
can be events, strategies, or policies as nodes of the graph and relationships as edges, where each 
concept is mathematically represented as a neutrosophic vector from the neutrosophic vector space. 
Every node, in its vector form, is represented by (x1, , Xn); x; € {0,1,1}, where 0 is off state, 1 is on 
state and I is the indeterminate state. 

Consider two nodes N, and N, of the NCM; their relationship is given by the edge e,,. Every 
weighted edge e,, is from {-1, 0, 1, 3, where 0 means no impact, positive value means increase in 
Na implies increase in Np, similarly decrease in N, implies decreases in Np. If egy takes a negative 
value like —1, it implies that an increase in Ng implies a decrease in Np, or similarly, a decrease in 
Nq implies an increase in N,. The edge weight is assigned a value I if the effect from N, to N, can 
not be determined. The edge weights of simple NCMs are from {-1, 0, 1, I}. The NCM’s adjacency 
matrix is denoted by N(E) = (e ap), Cay € {-1,0,1,]} . 
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NCMs are expert opinions constructed based on data obtained from the expert, where they 
identify the factors or concepts relevant to the domain and associated causal relationships in terms of 
numbers and indeterminacy or linguistic terms. 

Generally, the neutrosophic connection matrix is obtained directly from the digraph of the 
expert’s opinion. Here, we have introduced the concept of dealing with neutrosophic linguistic terms 
using the following logic. Algorithm 1 provides neutrosophic edge weights from the linguistic terms. 


Algorithm 1: Generate NCM Matrix 


Data: matrix « Empty Matrix; row «<- Empty List; n « No of nodes 
Result: NCM Matrix; 
i,j — 0 while i#ndo 
string + “” 
while 7 # n do 
term < input linguistic term 
if i= j then 
| element «+ 0 


else 

if term = ”"-VH” then 

| element <—- random(—5, —2.75) 
else if term = ”"-H” then 

| element «< random(—2.75, —2) 
else if term = ”-M” then 

| element <— random(—2, —1.5) 
else if term = ”-L” then 

| element «+ random(—1.5, —1) 
else if term = "-VL” then 


| element «+ random(—1,0) 
else if term = "NC” then 

| element — 0 
else if term = "+VH” then 

| element < random(2.75, 5) 
else if term = "+H” then 

| element < random(2, 2.75) 


else if term = "+M” then 
| element «<- random(1.5, 2) 
else if term = "+L” then 


| element < random(1, 1.5) 
else if term = "+VL” then 

| element <- random(0, 1) 
else if term = “NaN” then 

| element < I 


string < string + ele ment+' ,’ 
3t-5+1 
for « € temp do 
if «="J" then 
| | continue 
row ¢+- row.append(temp) 
#<i+1 
itO 
while i # n do 
matrix <~ matrix.row-insert(i, row{i]) 
i¢-i+1 
return matrix 


In the case of linguistic terms, these edges can be considered as negative very high (-VH), 
negative high (-H), negative medium (-M), negative low (-L), and negative very low (-VL). It implies 
a negative influence when the decrease of the influence of the node N, results in the increase of the 
influence of the node N, or the increase of N, results in the decrease of N,. Similarly, these edges 
can be considered as positive very high (+VH), positive high (+H), positive medium (+M), positive 
low (+L), and positive very low (+VL). It implies a positive influence when the decrease of N, results 
in the decrease of N, or the increase of N, results in the increase of Nj. If the linguistic term is no 
causality, 0 is used; if it is indeterminate, “NaN" (not a number) is used. These linguistic terms are 
converted to numerical values using simple assignments as follows: 

e If term = “-VH" then the element gets a random value from [-5, -2.75). 
e If term = “-H", then the element gets a random value from [-2.75, -2). 
e If term = “-M", then the element gets a random value from [-2, -1.5). 
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e If term = “-L", then the element gets a random value from [-1.5, -1). 
elf term = “-VL", then the element gets a random value from [-1, 0). 

e If term = “-NC" then the element gets 0. 

e If term = “+VH" then the element gets a random value from (2.75,5]. 


e If term = “+H", then the element gets a random value from (2, 2.75]. 
e If term = “+M", then the element gets a random value from (1.5,2]. 
e If term = “+L", then the element gets a random value from (1, 1.5]. 
elf term = “+VL", then the element gets a random value from (0, 1]. 


e If term = “NaN" then the element gets I. 

The values generated by the algorithm range from [-5, 5]; since the non-indeterminate edge 
weights of NCM are from [-1, 1], the edge weights are normalized before the NCM/neutrosophic 
digraph is constructed. 

The expert opinion obtained can be used to generate the neutrosophic adjacency matrix of the 
NCMs. Generation can be done by using linguistic terms. According to the linguistic terms, the 
neutrosophic matrix can be obtained using random values between the ranges. 


Example 1: Consider the graph of the NCM, using the seven concepts (or nodes or attributes) 
L,,L2,..,L7 and the expert opinion is obtained in terms of linguistic terms, using Algorithm 1, the 
weights for the edges are provided. 

For illustration, assuming that the connection from L, to L7 has a very high positive influence, 
and the connection between from L, to L3 is indeterminate, and the connection from L, to L, isa 
very high negative influence. The edge weights generated by the algorithm using the linguistic terms 
are tabulated in the Table 1: 


Table 1. The edge weights assigned by the algorithm. 


iy Ly i La i ie L, 
i 0 -25 I 15 0 3.5 0 
Ly 0 0 0 0 15 (OO 5 
be 0 0 0 0 5 2 0 
La 0 I 0 0 0 0 -45 
ie 0 0 0 0 0 0 eS] 
Le | 05 | O 0 0 0 0 I 
L 0 0 0 0 0 0 0 


The NCM will be given in Figure 1. Each edge is weighted and directed. The dashed lines are used to 
represent indeterminate edges. The weight of each edge is given in the graph. 


Figure 1. An illustration of the neutrosophic directed graph. 
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It is then normalized to weigh between [—1,1] and indeterminacy I. The connection matrix for the 


neutrosophic directed graph is given in Eq. (1). 


0 —05 J 03 0 0.7 0 
0 0 0 O 0.3 0 1 
0 0 0 O 1 04 O 
0 I 0 O 0 0 —0.9 
E normalized = 0 0 0 0 0 0 —0.2 (1) 
0.1 0 0 O 0 0 I 
0 0 0 O 0 0 0 


The edge weights can also be obtained from the expert. The edges of simple NCMs are from {-1, 0, 
1, B. 


Example 2: Consider a simple NCM given by an expert with seven concepts L,,L2,...,L7 as the nodes 
of the directed neutrosophic graph. The expert opinion is obtained in terms of edge weights. 


Figure 2. Neutrosophic directed graph for simple NCM. 


The edge weights are from {—1,0,1,/}, and the indeterminate edges are represented by dotted lines. 
The NCM’s connection matrix is denoted by N(E) =(e q,), where eg, € {0,1, —-1,/}. 


OO -1 7% £0 1 0 
0 0 0 0 1 0 1 
OF 0b. oD. ee 0 
oF 0 0 0 0 =1 

ve 0 0 0 0 0 0 -1 (2) 
Ld” i. 200 “OO? FF 
0 0 1 0 0 0 0 


The notion of state vector, dynamical system and its functioning are described. The neutrosophic 
state vector S = (sj,...,5)); 5; € {0,1,/}; where 0 indicates the off state, 1 is in the on state, and I 
indicates the indeterminate state. Let Egle: Likes, lala, ot rae be the NCM’s directed edges. 
Given the edges creating a cycle, the NCM becomes cyclic; otherwise, it is acyclic. Consider 1,13, 
Tsla. oe tea. te be cyclic, ifnode L, ison the influence will flow via the existing cycle and cause 
L, to be on again. This state of equilibrium of the dynamical system is called as a hidden pattern. 
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Various state vectors with different nodes in on state are considered to activate the system 
Consider the NCM with feedback given in Figure 2; its neutrosophic adjacency matrix is N(E) given 
in Eq. (2). 

The state vector S, = (1, 0, 0,, 0) where L, is in on state is considered. The data needs to be 
transformed by N(E), so we multiply S, by N(E). 

The state vector multiplied by the neutrosophic matrix N(E£) is given in Algorithm 3. In the 
given algorithm, S is a 1-dimensional matrix (row vector) and the neutrosophic adjacency matrix 
N(E) is a 2-dimensional matrix denoted by B. The resultant vector Res of the multiplication of S x 
N(E) is returned. 

An example illustrates this: Consider the graph in Figure 2 and its related connection matrix. 
Take the state vector S; = (1000000). The state vector S,; is multiplied with the neutrosophic 
adjacency matrix N(E). 


6 =i f ££ @ 2 6 
00 0 0 1 0 1 
0 0 00 1 1 «0 
= 0 F 0 0 © 0 -1 
Sx N(E)=(1000000)X], 9 9 9 9 0 <4 (3) 
11 00007 
> 0 1 0 6 0 Oo 
=(0 —1/1010) (4) 


Algorithm 2: Multiply 


Data: 1D Matrix S and B + 2D Neutrosophic Adjacency Matrix N(E) 
Result: Resultant Matrix 
Res — 1D Matrix 
r< S.size 
R<€ B.size 
i,j +0 
while i+ r do 
while j 4 R do 
Res{i] — Res{i] + (S{i] *« B[j][i)) 
j¢jtl 
i<-it+l 
return Res 


After obtaining the resultant vector Res, it must undergo the threshold and update operations. By the 
definition of NCM, the threshold and update operation is denoted by the © symbol. 

It is important to note here that working with I (indeterminate) needs to be done carefully. As 
by the definition of I 


Ixl=P=1 (5) 


Any power of I gets mapped to I, as shown in Algorithm 3. 


Algorithm 3: UpdateIPower 
Data: V + state vector 
Result: V state vector after updating /’s 
i+O 
forie€V & j€bdo 
| ifi=J? then ic I 
return V 


Similarly, during the updating and threshold operation, any constant into I is also mapped into I, 
as shown in Eq. (6) 
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nxl=I (6) 


The resultant vector from the multiplication of the state vector with N(E) is thresholded and 
updated. 
Let X =S ,N(E) =(s 1,5 2,, 8») is thresholded by replacing s; accordingly to the Eq. (7). 


1 if Si >t 
s, = 40 if s;<t (t is a suitable positive integer) (7) 
I if s;not an integer 


The resultant X is updated to ensure that the state considered on in the initial state vector S, is on 
in the resulting vector. Here, it is updated to make the concept N , as 1 in the resulting vector. 
The algorithm for thresholding and updating is given in Algorithm 4. 


Algorithm 4: ThresholdandUpdate 


Data: X < resultant vector and t threshold value, Initial On state 
Result: Thresholded and updated resultant vector 
X < updateIPower(X) 
len + X.size 
i+O0 
while i 4 len do 
temp-_expr < X/i| 
if temp_ecpr = J then 
| temp_expr + 0) 
if temp_expr > t then 
| X[ij<t 
else if temp_expr = 0 then 
if X{i] 40 then 
| X[iJ<cJ 
else 
| X[i] <0 
4<-t+1 
X(state-—1)<+ 1 
return X 


Here, the Algorithm 4 illustrates where only one node is taken in the on state in the initial state vector. 
The thresholding and updating operation is mathematically denoted by ©. For example, 
consider the resultant vector of S, x N(E), the thresholding and updating result in S). 


S, x N(E) = (0 —111010) 4 (1011010) =S, (8) 


Here, during the threshold operation, —1 is made 0, and during the update operation, the very 
first state is made on again. 

Generally, in any NCM, more than one concept/ node can be considered in the on state. This will 
deal with the combined effect of both states being on. The proposed model can handle more than one 
node in the on state, and it can work to analyse the effect of a combination of various nodes. The 
multiplication of resultant vectors with the neutrosophic adjacency matrix N(E£) will continue until 
the resultant vector yields a fixed point or limit cycle. 

If the NCM settles to a neutrosophic state vector repeating in the form 


$198,405 «48; 4S, 


Then the dynamic system’s equilibrium is called NCM’s limit cycle. Suppose, S ,N(E) & S, 
(where © denotes the resultant vector of S,;N(E) which is thresholded and updated) and for S ,N(E) 
we repeat the same procedure until we attain the fixed point / limit cycle. 
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Two vectors are compared in Algorithm 5. 


Algorithm 5: Compare 


Data: a,b - 1D vectors 
Result: Boolean value 
i,j <0 
forica & jE bdo 
if i+ j then 
| return false 
return true 


A fixed point or limit cycle is attained when both the vectors under comparison are the same. 
For example, if we compare S, with Sj, we can see that S, # S. 


(1000000) #(10/1010)=S, (9) 
Since the fixed point nor limit cycle is reached, i.e., 5; # Sz, the process is continued. 
S.X N(E) = (11111141 -14+D°(1111110) =S3 (10) 
Since the fixed point still needs to be reached, i.e., S3 # Sz, the process is continued. 


Sy X N(E) =(1+1 27 21.1 22147 -14+2/ 
% (1111110)=S, (11) 


Here S3 = S,. The fixed point has been reached. 


The Algorithm 6 is used to determine if a limit cycle is reached. It compares with the previous 
resultant vectors using the previously described in compare Algorithm 5 to check if the limit cycle is 
reached. In case it is reached, it returns a true or a false. 


Algorithm 6: Check_cycle 
Data: b < list of vectors 
Result: Boolean value 
n < b.size 
i+O 
while i#n—1do 

jecitli 

while j #4 do 
vl + Di] 
v2 + b{j] 
if compare(vl,v2) = true then 

| return true 

jejt+l 

i¢it+l 

return false 


To illustrate an example of limit cycle, consider the state vector P; = (0100000), and the 
connection matrix N(E) as given in Eq. (2), that is the related adjacency matrix for the NCM given in 
Figure 2. 


P, x N(E) = (0000101) © (0100101) =P, 
P, X N(E) = (0010100) 6(0110100) =P, 
P; x N(E) = (0000210) 6(0100110) =P, 
P, x N(E) = (1100107) 9 (1100101 =P; 
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P; x N(E) = (0-12*11110) 
6 (0111110) =P, 
P, x N(E) =(11+1001+111-1) 
6 (1100110) =P, 
P, x N(E) = (7I1-111111?) 
6 (1111111 =P, 
Pg X N(E) = (11/2241 1 1+12*I1-1) 
6 (1111110) =P, 
PyX N(E) = (2 +1-11114114+1 2-1 
6 (1111110) =P, 
Py X NE) = (1122711 4+12*11-1) 
6 (1111110) =Py, = Py (12) 


Since P,, = Py, the iteration is stopped since a limit cycle has been achieved, enabling the 
determination of the hidden pattern. The limit cycle is as follows: Pg gives Py, Py gives Pio, Pio 
gives P,,; thatis same as Po. 


The process of multiplication of the resultant vector with the matrix N(E) is continued until a 
limit cycle / fixed point is reached. 

The Algorithm 7 takes an adjacency/connection matrix and state vector (an integer as input 
denoting the state to be activated) and an integer denoting the threshold values as parameters. 


Algorithm 7: Iteration 
Data: an Adjacency matrix EF, state < state to be activated, t threshold value 
Result: list of vectors 
len + E.size 
cl < list{len] filled with 0 
cl{state —1] + 1 
start + Matriz(cl) 
flag < false 
vectors < Empty List 
while flag = false do 
y + multiply(start, E) 
y + thresholdAndUpdate(y, t) 
vectors < vectors.append(y) 
start + y 
flag + check_cycle(vectors) 
return vectors 


The algorithm 7 continues multiplying the resultant vector with the N(E) until a fixed point ora 
limit cycle is reached. It is dependent on several previously described algorithms. Every concept must 
be made in the active state to capture the hidden pattern and understand the effect of the concept or 
node on others to analyse the NCM thoroughly. Only in NCMs can we signify that the clout of anode 
on different nodes can be indeterminate, and this vision needs to be revised in the case of FCMs. 

For the given example, the results have been tabulated in Table 2 for various state vectors. 
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Table 2. Determining the Hidden Pattern. 


Input state vector __ limit cycle / fixed point Resultant vector 
(1000000) fixed point (1171/10) 
(0100000) limit cycle (1111110) 

(1111110) 
(1111110) 
(1111110) 

(0010000) fixed point (111110) 
(0001000) fixed point (01011 00) 
(0000100) fixed point (0000100) 
(0000010) fixed point (1111110) 
(0000001) fixed point (711111) 


It is seen that only the active state of concept Lz results in a limit cycle—the rest results in a fixed 
point. 

The conclusions that can be drawn from the equilibrium state of the dynamic system are as 
follows: When node L, is active, all nodes are either indeterminate or on state. When L, isin the on 
state, it results in a limited cycle, which affects all nodes other than L,. Similarly, when nodes 
L3,L4,L¢ or Lz alone are in the on state, other nodes are either indeterminate or on state. Whereas 
when node L, is in on state, no other node is affected; all of them remain in the off state. 


3. Description of the modelling package 


The overall flow of the modelling software is as given in Figure 3. The first module is for the 
input module, which can either be linguistic term-based or edge-weight-based. 


Causality in Linguistic 
Terms 
« Assign Edge weights 
according to terms 


Creating the 


< Generate < 
Input No of Neutrosophic P Simulate the 
P : Connection : Results 
Graph using the : dynamic system 
Nodes : Matrix 
edge Weights 
ieee "State vector 
(OR) | * S2= S1*N(E) 
? \ 17,oos00 
r’ ~~ 1000000 * Threshold and update 
Direct Edge Weights ieee a , ave a operation 
from Expert opinion ¥, — 0000 ° * Limit Cycle/Fixed point 


Figure 3. Various modules of NCM modelling software. 


1. Input module: There are two methods in which the model can be created using expert opinion. 
They can enter by either method as described below: 

- Linguistic terms: The user can enter the linguistic terms to describe the relationship between 
two concepts in NCM as very negative or positive. According to the previously discussed 
algorithm, the edge weights are assigned based on the linguistic term. 

- Edge weights: The user can instead directly enter the edge weight to denote the causality 
between two nodes. 

2. Neutrosophic digraph: The neutrosophic digraph is generated and visualised using the edge 
weights obtained from the user. 
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3. Connection matrix: The related connection matrix N(E£) is obtained from the neutrosophic 
digraph of the NCM. 
4. Dynamical system: Using various state vectors, the dynamical system is simulated to analyse 
the effect of the on state of the various nodes. 
5. Results: The effect of the on state of various nodes and combination of various states is 
consolidated. 
The NCMpy python package performs a simulation of the dynamic system —a detailed description 
of the package is given in the next subsection. 


3.1 Description of the NCMPy Python Package 


The flowchart of the NCMpy package is given in the following Figure 4. 


I ompa: 
ee State Vector updateIPower () S zat? 


| | | 7 


Resultant Resultant Fixed Point/ 
= ae —yYes—> aoe 
Vector > thresholdAndUpdate () Vector check_cycle() Limit Cycle 


| . | 


Figure 4. Flowchart of the NCM package. 


iteration () ——> multiply () — 


The SymPy library has been used to handle indeterminate values in modelling the NCMs in this 
proposed package. In SymPy package, the names E, I, O, S, N, and Q collide with names already 
defined in the package. Hence, J can not be used to represent indeterminacy, so J is used instead of 
I. Throughout the coding snippets, J is used instead of I. 

The various functions/modules are described here. 


1. compare(x, y): This function takes two vectors as input and checks for their equality; it is 
dependent on Algorithm 6. The code snippet is as follows: 


def compare(x, y): 
s = True 
for k in range(len(x)): 
if x[k] != y[k]: 
return False 
return True 


Assume that compare function takes S; = (J]II110) and Sg =(IIII1I10) as x and y. In 
that case, compare(x, y) will return True, in case it was Ss; = (1/11/11) and S,=(111110) 
under consideration, then it would return False. 


2. check_cycle(b): This function is dependent on Algorithm 6; it is used to check if the resultant 
vector is a fixed point or limit cycle. 
def check_cycle(b) : 
for i in range(len(b) - 1) : 
for j in range(i + 1, len(b)) : 


v1 = b[i] 
v2 = b[j] 
if(compare(v1 , v2) == True) : 


return True 
return False 
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The fixed point is achieved when the resultant vector is the same as the previous resultant vector, 
that is P;N(E) © P;. The limit cycle is achieved when the recently calculated state vector is the same 
as any one of the previously calculated resultant vector, which results in a cycle. 


P:N(E) & Pita; PitiN(E) © Pita; PeN(E) & Pi 


In Example 2, considering the active state P; = (0100000) results in a limit cycle. 


3. updatelPower(x): This function takes a 1D vector as input and converts the indeterminate 
quadratic polynomial to a linear polynomial. This function is based on Algorithm 3. This function is 
used by the threshold and update function. 
def updateIPower(x): 

for i in range((np.shape(x))[1]): 
if x[i].find(3**2): 
x[i] = x[i].xreplace({J**2: J}) 
return x 


Consider the resultant vector; 
R, = CU **2JJ +110] ** 22/) 
For a sample scenario. The updateIPower(R_1) will change this vector P, into 


R, = (JJ +110] 2) 


4. thresholdAndUpdate(X,threshold_value,state) : This function takes a 1D vector and a 
threshold value as a parameter and updates each vector value according to the defined thresholding 
operation. Also, the updation operation must see to it that the node which was on in the initial state 
is on in the next state and so on in the resultant state also; if not, it is set to 1 again. 

def thresholdAndUpdate(X , threshold_value, state) : 
X=updateIPower(X) 
for i in range((np.shape(X)[1])): 
if(X[i].find(-1)): 
X[i]=e 
temp_expr=X[i].subs(J,@) 
if (temp_expr>=threshold_value): 
X[iJ]=X[i].subs(X[i],threshold_value) 
elif(temp_expr==0): 
if(X[i].find(J)): 
X[i]=3 
else: 
X[i]=X[i].subs(X[i],@) 
if isinstance(state, list): 


activated_states = [i for i, x in enumerate(state) if x == 1] 
for s in activated_states: 
X[s] =1 
else: 
X[state - 1] =1 
return X 


Xi XN(E)= (F417 +1777 412«11-1)9 (1111110) =X, 


5. iteration(E,state,threshold_value) : This function takes an adjacency/connection matrix, an 
integer denoting the state to be activated and an integer denoting the threshold values as parameters. 
This function is based on the Algorithm 7 . 
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def iteration(E , state, threshold_value = 1) : 
if isinstance(state, list): 
start = sym.Matrix(state) 
else: 
cl = np.zeros((np.shape(E)[1])) 
ci[state - 1] =1 
start = sym.Matrix(c1) 


flag = False 

start = start.T 

vectors = [] 

while flag == False : 
y = multiply(start , E) 
# Performing the thresholding operation on output vector y 
y = thresholdAndUpdate(y , threshold_value, state) 
vectors. append(y) 
# Updating start vector to start with new state vector 
start = y 
# Checking for cycle among state vectors 
flag = check_cycle(vectors) 

return vectors 


Consider Example 2, where the working out is done with active state P; = (0100000); it 
iterates until the limit cycle is achieved. 


3.2 Modelling NCMs for sample case study 


The opening page for the modelling software is given in Figure 5. Here, the user can select the 
option of working with linguistic terms or directly entering the neutrosophic edge weights given by 
the expert as shown in Figure 6. 


WE HELP YOU TO 


Model Your 
NCMs 


Figure 5. Homepage. 


How would you like to give input? 


Figure 6. Option selection. 


3.3 NCMs using Linguistic Terms 


The neutrosophic linguistic terms are obtained from the user, as shown in Figure 7. 
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t Linguistic Terms 


npu 
BE eer tne number ot nodes Add Nodes 


+VH +H -HO +H +M +H 
+L NON NGN +H O -L NaN 
-L-H +VH +M-LO +VH 


-H -VH +M +H NaN +VHO 


Figure 7. Sample input page for the linguistic terms related to the neutrosophic directed graph. 


The NCMpy package runs through the Alogrithm GenerateNCM 2 and creates the necessary 
edge weights and normalises them as shown in Table 3, providing the neutrosophic bigraphs. The 


neutrosophic-directed graph of the NCM is given in Figure 8. 


Table 3. The edge weights assigned. 


Ls Ly i i; te 
im 0 4.88 I 3.58 1.91 
i 4.22 0 I 2.13 I 
ls I I 0 -4.77 I 
La 2.98 158 -2.41 0 2.41 
be 133 I I 4.2 0 
i. | 120.436.) 4.03.) d71 | --240 
L, | -2.03 | -4.77 | 1.34 | 3.89 I 


The resultant neutrosophic connection matrix of the graph is 


[ 

[@, 0.98, J, 0.72, 0.38, -@.48, -0.56], 
[0.84, @, J, 0.43, J, -0.55, -@.94], 

[J, J, @, -0.95, J, 0.65, 0.39], 

[0.60, 0.32, -0.48, 0, 0.48, 0.36, 0.62], 
[0.24, J, J, 0.84, 0, -0.42, J], 

[-0.24, -0.87, 0.81, 0.34, -0.48, @, 0.57], 
[-0.41, -0.95, 0.27, 0.78, J, 0.89, @] 

] 
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0.72 


0.39 


0.24 


Figure 8. Neutrosophic directed graph. 


To show the plotting capacity of the visualizing module, we have taken a matrix with all connections 
for the sample. The threshold value is obtained from the user. According to the threshold value set 
by the user, the thresholding and updating of state vectors are done. 


Figure 9. Insert the threshold value and state vector. 
The threshold value of 1 was taken here, and the following results were obtained for various state 


vectors. The working out for the first state vector S$, = (1000000), where the concept C, is on is 
shown here: 


Kandasamy et al., NCMPy: A Modelling Software for Neutrosophic Cognitive Maps based on Python Package 


Neutrosophic Systems with Applications, Vol. 13, 2024 17 
An International Journal on Informatics, Decision Science, Intelligent Systems Applications 


FOR ACTIVE STATE 1 

Matrix([[@, 0.98, 1.0*J, 0.72, 0.38, -@.48, -0.56]]) 

Matrix([[@, @.98, 1.0*J, 0.72, 0.38, -0@.48, -0.56]]) 

Matrix([[@, 0.98, 1.0*J, 0.72, 0.38, -@.48, -0.56]]) 

Matrix([[@, @, 1.0*J, 0.72, 0.38, -@.48, -@.56]]) 

Matrix([[@, @, J, 0.72, 0.38, -@.48, -@.56]]) 

Matrix([[@, @, J, 0, 0.38, -@.48, -@.56]]) 

Matrix([[@, @, J, @, @, -@.48, -0.56]]) 

Matrix([[@, @, J, @, @, @, -0.56]]) 

After Thresholding and Updating 

Matrix([[1, ®, J, ®, @, @, @]]) 

Matrix([[3**2, J**2 + 0.98, J, @.72 - @.95*J, J**2 + 0.38, @.65*] - 0.48, @.39*] - @.56]]) 
Matrix([[J, J + 0.98, J, @.72 - @.95*J, J + 0.38, @.65*J - 0.48, 0.39*J - 0.56]]) 
Matrix([[J, J] + 0.98, J, 0.72 - @.95*J, J + 0.38, @.65*J - 0.48, @.39*J - @.56]]) 
Matrix([[J, @, J, 0.72 - @.95*J, J + 0.38, @.65*J - 0.48, @.39*J - 0.56]]) 
Matrix([[J, @, J, @.72 - @.95*J, J + @.38, 0.65*J - @.48, @.39*) - @.56]]) 
Matrix([[J, @, J, @, J + 0.38, @.65*J - 0.48, 0.39*J - @.56]]) 

Matrix([[J, 0, J, @, @, @.65*J - 0.48, 0.39*J - @.56]]) 

Matrix([[J, ®, J, @, @, 0, @.39*J - @.56]]) 

After Thresholding and Updating 

Matrix([[1, ®, J, @, @, @, @]]) 

[Matrix([[1, @, J, 0, @, @, Q}]), Matrix([[1, @, J, 0, 9, 9, 0}]1)] 


Similarly, the working out for each and every state vector is carried out. For state vector By = 
(0 100000), where the concept C, is on, the resultant vectors will be: 


FOR ACTIVE STATE 2 

Matrix([[0.84, 0, 1.0*J, 0.43, 1.0*J, -0.55, -@.94]]) 

After Thresholding and Updating 

Matrix([[®, 1, J, @, J, @, @]]) 

Matrix([[J**2 + 0.24*J + 0.84, 2*]**2, J**2 + J, 0.43 - @.11*J, J**2 + J, @.23*J - 0.55, J**2 + @.39*J - @.94]]) 
After Thresholding and Updating 

Matrix([[@, 1, J, @, J, @, @]]) 

Resultant vector(s) [Matrix([[@, 1, J, @, J, 0, @]]), Matrix([[@, 1, J, 0, J, @, @]])] 


For state vector G,; = (0010000), where the concept C3 is on, the resultant vectors will be: 


FOR ACTIVE STATE 3 

Matrix([[1.0*J, 1.0*J, @, -0.95, 1.0*J, 0.65, @.39]]) 

After Thresholding and Updating 

Matrix([[J, J, 1, 0, J, @, @]]) 

Matrix([[2.1*J, J**2 + 2.0*J, 3*J**2, 2.0*J - 0.95, J**2 + 1.4*J, 0.65 - 1.4*J, J**2 - 1.5*J + 0.39]]) 
After Thresholding and Updating 

Matrix([[J, J, 1, @, J, @, @]]) 

Resultant vector(s) [Matrix([[J, J, 1, ®, J, 0, 0]]), Matrix([[J, J, 1, ®, J, 0, @]])] 


For state vector X,; = (0001000), where the concept C, is on, the resultant vectors will be: 


FOR ACTIVE STATE 4 

Matrix([[@.60, @.32, -0.48, 0, 0.48, 0.36, 0.62]]) 

After Thresholding and Updating 

Matrix([[@, ®, @, 1, ®, @, @]]) 

Matrix([[@.60, 0.32, -@.48, 0, 0.48, 0.36, @.62]]) 

After Thresholding and Updating 

Matrix([[@, @, @, 1, @, @, 0]]) 

Resultant vector(s) [Matrix([[@, ®, @, 1, ®, 0, @]]), Matrix([[@, 9, 0, 1, 0, @, @]])] 


For state vector Y, = (0000100), where the concept C; is on, the resultant vectors will be: 


FOR ACTIVE STATE 5 

Matrix([[@.24, 1.0*J, 1.0*J, 0.84, @, -@.42, 1.0*J]]) 

After Thresholding and Updating 

Matrix([[@, J, J, @, 1, @, J]]) 

Matrix([[J**2 + 0.43*] + 0.24, J**2 + @.05*J, J**2 + 1.3*J, 0.26*J + 0.84, 3*J**2, @.99*J - 0.42, @.45*J]]) 
After Thresholding and Updating 

Matrix([[@, J, J, ®, 1, @, J]]) 

Resultant vector(s) [Matrix([[@, J, J, @, 1, ®, J]]), Matrix([[@, J, J, @, 1, @, J]])] 


For state vector Z, = (0000010), where the concept C, is on, the resultant vectors will be: 


FOR ACTIVE STATE 6 

Matrix([[-@.24, -0.87, 0.81, 0.34, -0.48, @, 0.57]]) 

After Thresholding and Updating 

Matrix([[@, @, @, @, @, 1, @]]) 

Matrix([[-@.24, -0.87, 0.81, 0.34, -0.48, @, 0.57]]) 

After Thresholding and Updating 

Matrix([[@, @, @, @, @, 1, @]]) 

Resultant vector(s) [Matrix([[@, ®, @, 0, @, 1, 0]]), Matrix([[®, 0, ®, @, 0, 1, @]])] 
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FOR ACTIVE STATE 7 
Matrix([[-@.41, -@.95, 0.27, @.78, 1.0*J, 0.89, @]]) 
After Thresholding and Updating 
Matrix([[@, @, @, @, J, @, 1]]) 
Matrix([[@.24*] - @.41, J**2 - 0.95, J**2 + 0.27, @.84*J + 0.78, J, 0.89 - 0.42*J, J**2]]) 
After Thresholding and Updating 
Matrix([[@, 0, ®, ®, J, ®, 1]]) 
Resultant vector(s) [Matrix([[@, 9, ®, ®, J, @, 1]]), Matrix([[@, 0, ®, @, J, @, 1]])] 
For state vector A; = (000000 1), where the concept C; is on, the resultant vectors will be: 
The result vector in each case has been shown. Results regarding the resultant vectors that can 
be discussed 
1. Maximum Influence: The nodes C,, C3; and C; are the most influential since they do 
affect many other nodes and turn to an indeterminate state. 
2. Least Influential nodes: The nodes C, and C, are the least influential since they do not 


affect any other node than itself. 


3.4 NCMs using Edge Weights 


The edge weights are obtained from the expert, as shown in Figure 10. 


Input Neutrosophic Matrix 


ibor Political-Leaders Goog 


Figure 10. Input screen. 


Using the edge weights obtained, the NCM is created. It is the visualization of the same represented 
as a NCMs digraph as shown in Figure 11. 


Figure 11. Visualization of the neutrosophic graph. 


The related neutrosophic connection matrix N(E£) is given in Eq. (2). The obtained matrix is 


(9, -1, "J", 1, @, 1, Q], 
(9, 9, , @, 1, @, 1], 
(0, ®, @, @, 1, 1, 9], 
[@, ‘J*, 0, 0, @, ®, -1], 
(9, 0, @, @, ®, 0, -1], 
[1, 1, ®, @, @, 0, 'J'], 
[9, @, 1, @, @, @, @] 
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The threshold value of 1 was taken here, and the following results were obtained for various 
state vectors. The working out for the first state vector S, = (1000000), where the concept C, is 
on is shown here: 


FOR ACTIVE STATE 1 


Matrix([[@, -1.ee@eeeeeeeeeeee, 1.0*3, 1.e@e@@@eeGe0R00, 0, 1.eee0eeeee0G000, 0]]) 


After Thresholding and Updating Matrix([[1, @, J, 1, 0, 1, @]]) 

Matrix([[1, J, J, 1, J, J +1, J - 1)]) 

After Thresholding and Updating Matrix([[1, J, J, 1, J, 1, @]]) 

Matrix([[1, J, J, 1, 2*J, J +1, 3 - 1]]) 

After Thresholding and Updating Matrix([[1, J, J, 1, J, 1, @]]) 

Resultant vector(s) [Matrix([[1, ®, J, 1, ®, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]])] 
The resultant vectors for all active states are given in the snippet below: 

FOR ACTIVE STATE 2 

Resultant vector(s) [Matrix([[@, 1, @, @, 1, @, 1]]), Matrix([[@, 1, 1, @, 1, ®, @]]), Matrix([[®, 1, @, 

®, 1, 1, @]]), Matrix([[1, 1, ®, @, 1, @, J]]), Matrix([[@, 1, J, 1, 1, 1, @]]), Matrix([[1, 1, ®, @, 1, 

J, @]]), Matrix([[J, 1, J, 1, 1, 1, J]]), Matrix([[1, 1, J, J, 1, J, @]]), Matrix([[J, 1, J, 1, 1, 1, 

@]]), Matrix([[1, 1, J, J, 1, J, @]])] 

FOR ACTIVE STATE 3 

Resultant vector(s) [Matrix([[@, @, 1, @, 1, 1, @]]), Matrix([[1, 1, 1, ®, 1, 1, @]]), Matrix([[1, 9, 1, 

1, 1, 1, J]]), Matrix([[1, J, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, @]])] 

FOR ACTIVE STATE 4 

Resultant vector(s) [Matrix([[@, J, @, 1, 0, @, @]]), Matrix([[@, J, @, 1, J, @, @]]), Matrix([[@, J, @, 

1, J, @, @]])] 

FOR ACTIVE STATE 5 

Resultant vector(s) [Matrix([[@, @, @, @, 1, @, @]]), Matrix([[@, @, @, @, 1, 9, @]])] 

FOR ACTIVE STATE 6 

Resultant vector(s) [Matrix([[1, 1, @, @, @, 1, J]]), Matrix([[1, @, J, 1, 1, 1, 1]]), Matrix([[1, J, 1, 

1, J, 1, @]]), Matrix([[1, J, J, 1, 1, 1, @]]), Matrix([[1, J, J, 1, J, 1, 0]]), Matrix([[1, J, J, 1, J, 

1, @]])] 

FOR ACTIVE STATE 7 

Resultant vector(s) [Matrix([[@, @, 1, @, ®, @, 1]]), Matrix([[@, @, 1, @, 1, 1, 1]]), Matrix([[1, 1, 1, 

@, 1, 1, 1]]), Matrix([[1, ©, 1, 1, 1, 1, 1]]), Matrix([[1, J, 1, 1, 1, 1, 1]]), Matrix([[1, J, 1, 1, 1, 


1, 1)))] 


The combination of two concepts being in the on state is also handled by the package done there. 
Similarly, several other combinations can be done, according to the user’s choice. A sample of the 
combinations of various concepts begins in on state, and the working is given in the snippet. 


Activating state: [i; 1; 0; 0; 0, 0, e] 

(Matrix([[1, 1, J, 1, 1, 1, 1]]), Matrix([[1, 1, 1, 1, 1, 1, @]]), Matrix([[1, 1, J, 1, 1, 1, @]]), 
Matrix([[1, 1, J, 1, 1, 1, ©]])] 

Activating state: [1, 0, 1, @, 0, 8, @] 

[Matrix([[1, @, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, @]]), Matrix([[1, J, 1, 1, 1, 1, @]])] 
Activating state: [1, @, @, 1, 9, @, @] 

[Matrix([[1, @, J, 1, @, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]]), Matrix([[1, J, J, 1, J, 1, ]])] 
Activating state: [1, 0, @, @, 1, 98, @] 

[Matrix([[1, @, J, 1, 1, 1, @]]), Matrix([[1, J, J, 1, 1, 1, @]]), Matrix([[1, J, J, 1, 1, 1, @]])] 
Activating state: [1, 0, 0, @, @, 1, @] 

[Matrix([[1, @, J, 1, @, 1, J]]), Matrix([[1, J, J, 1, J, 1, @]]), Matrix([[1, J, J, 1, J, 1, @]])] 
Activating state: [1, @, @, 0, 0, @, 1] 

[Matrix([[1, ®, 1, 1, 6, 1, 1)]), Matrix([[1, J, 1, 1, 1, 1, 1))), Matrix([[1, J, 1, 1, 1, 1, 1)))] 
Activating state: [@, 1, 1, 0, @, @, @] 

[Matrix([[@, 1, 1, @, 1, 1, 1]]), Matrix([[1, 1, 1, ®, 1, 1, J]]), Matrix([[1, 1, 1, 1, 1, 1, J]]), 


Matrix([[15 1, 1, 2) 25 25. @]])s Matrax([[1, 15. 25 25 25.2, 01])] 


The complete NCMPy python package will be made available online. 


4. Conclusions and discussions 


As the world moves towards no code or less coding paradigms, a Python package for NCM will 
aid and help mathematicians, social scientists, economists, strategists, and other policymakers 
analyze real-world problems without worrying about the mathematical background or 
computational complexities of NCM. A visualization tool and a Python package to help in the 
working of NCM were presented in this paper. The NCMpy package and modelling software provide 
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the functions essential to studying problems involving indeterminacy, which can be done using 
NCMs. 

This package and modelling tool are open-source, written in Python, straightforward to 
implement, and provide the required functionality for handling models with indeterminacy. 

This tool implementation is a collaboration with the founding and leading experts in the field of 
NCMs. This tool will facilitate research and enable new researchers and scientists to apply NCMs to 
their projects that involve indeterminacy. We plan to update our library and constantly welcome all 
scientific community contributions. 
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